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■ 1 Introduction 

^ ■ The thermal conductivity of insulating (dielectric) crystals is computed almost ex- 

' clusively on the basis of the phonon Boltzmann equation. We refer to jl] for a 

^ . discussion more complete than possible in this contribution. On the microscopic 



X 



level the starting point is the Born-Oppenheimer approximation (see Plj for a mod- 



O 

O I ern version), which provides an effective Hamiltonian for the slow motion of the 
^ ' nuclei. Since their deviation from the equilibrium position is small, one is led to 
^ ' a wave equation with a weak nonlinearity. As already emphasized by R. Peierls in 
p |. his seminal work [Hj, physically it is of importance to retain the structure resulting 
fH ! from the atomic lattice, which forces the discrete wave equation. 

! On the other hand, continuum wave equations with weak nonlinearity appear in 

the description of the waves in the upper ocean and in many other fields. This topic 
is referred to as weak turbulence. Again the theoretical treatment of such equations 
is based mostly on the phonon Boltzmann eqation, see e.g. 0. In these applications 
one considers scales which are much larger than the atomistic scale, hence quantum 
effects are negligible. For dielectric crystals, on the other side, quantum effects are of 
importance at low temperatures. We refer to [1] and discuss here only the classical 
discrete wave equation with a small nonlinearity. 

If one considers crystals with a single nucleus per unit cell, then the displacement 
field is a 3- vector field over the crystal lattice F. The nonlinearity results from the 
weakly non-quadratic interaction potentials between the nuclei. As we will see, the 
microscopic mechanism responsible for the validity of the Boltzmann equation can 
be understood already in case the displacement field is declared to be scalar, the 
nonlinearity to be due to an on-site potential, and the lattice F = Z'^. This is the 
model I will discuss in my notes. 

As the title indicates there is no complete proof available for the validity of 
the phonon Boltzmann equation. The plan is to explain the kinetic scaling and to 
restate our conjecture in terms of the asymptotics of certain Feynman diagrams. 
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2 Microscopic model 

We consider the simple cubic crystal 1? . The displacement field is denoted by 

g^eM, a;eZ^ (2.1) 

with the canonically conjugate momenta 

p^eM, a;eZ^ (2.2) 

We use units in which the mass of the nuclei is m = 1. The particles interact har- 
monically and are subject to an on-site potential, which is divided into a quadratic 
part and a non-quadratic correction. Thus the Hamiltonian of the system reads 

^=^E (^''+^0?')+^ E «(x-y)?.g,+ Enfe) = ^o+E^(?-)- (2-3) 

The coupling constants have the properties 

a{x) = a[-x) , (2.4) 

\a{x)\ < aoe-^^""^ (2.5) 

for suitable ao, 7 > 0, and 

E«(x)=0, (2.6) 

because of the invariance of the interaction between the nuclei under the translation 

For the anharmonic on-site potential we set 

V{u) = Vl\xu^ + eiX^/lSuDu"^ , ueR. (2.7) 

s is the dimensionless scale parameter, eventually £ — > 0. The quartic piece is added 
so to make sure that H > 0. In the hmit £ — > its contribution will vanish and for 
simplicity of notation we will omit it from the outset. Then the equations of motion 
are 



x)qy{t) - ulq^{t) - y/EXqa^itf , x G . (2.8) 

We will consider only finite energy solutions. In particular, it is assumed that 
\Px\ 0, \qx\ — > sufficiently fast as |a;| 00. In fact, later on there will be the 
need to impose random initial data, which again are assumed to be supported on 



d , , 
Jt 



Px{t) , 
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finite energy configurations. In the kinetic limit the average energy will diverge as 
It is convenient to work in Fourier space. For / : — > R we define 



£-3 



f{k) = e"''"'"/. , (2.9) 



e = [-1 with inverse 



= / dke'^-'-^fik) , (2.10) 

dk the 3- dimensional Lebesgue measure. The dispersion relation for the harmonic 
part Ho is 

u;{k) ^ {ujl + a{k)Y^^ >u;o>0, (2.11) 

since a{k) > for A; 7^ because of the mechanical stabihty of the harmonic lattice 
with vanishing on-site potential. 

In Fourier space the equations of motion read 

d 

— g(A;, t) = p{k, t) , 
^^p(k,t) ^ -u;(kfq(k,t) 

-y/eX / dkidk2S{k - ki - k2)q{ki, t)q{k2, t) (2.12) 

with k eT^. Here 5 is the 5-function on the unit torus, to say, S{k') carries a point 
mass whenever k' El?. 

It will be convenient to concatenate and into a single complex- valued field. 
We set 

a{k) = l^{^q{k)+i^p{k)) (2.13) 

with the inverse 

m = ^-^Hk) + a{-kr) , m = ^^V^( - a(^) + a{-kr) . (2.14) 

To have a concise notation, we introduce 

a{k, +) = a{ky , a{k, -) = a{k) . (2.15) 

Then the a-field evolves as 
d 



{k,a,t) — iauj{k)a{k,a,t) + ia^/eX / dkidk2 

{8uj{k)u){ki)u){k2))~^^'^S{—ak + aiki + a2k2)a{ki, ai, t)a{k2, (72, t) . (2.16) 
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3 Kinetic limit and Boltzmann equation 



The kinetic limit deals with a special class of initial probability measures. Their 
displacement field has a support of linear size and average energy of order 
e~^. More specifically, these probability measures have the property of being locally 
Gaussian and almost stationary under the dynamics. Because of the assumed slow 
variation in space the covariance of such probability measures changes only slowly, 
i.e. on the scale e~^, in time. 

Let us assume then that the initial data for ()2.1(jj) are random and specified by 
a Gaussian probability measure on phase space. It is assumed to have mean 

{aik,a))^ = (3.1) 

and for the covariance we set 

{a{k,a)a{k',a))'^ = 0, (3.2) 

W'{y, k)=e^ [ dr]e^^y''{a{k - er]/2, +)a{k + er]/2, -))f , (3.3) 

y G {s'LY, which defines the Wigner function rescaled to the lattice (eZ)^. Local 
stationarity is ensured by the condition 

hm W {[r \,, k) = W°{r,k), (3.4) 
where [rj^ denotes integer part modulo e. Note that is normalized as 

V / dkW'{y,k)= [ dk{a{k,+)aik,-))^ . (3.5) 



The condition that the limit in ()3.4|) exists thus implies that the average phonon 
number increases as e~^, equivalently the average total energy increases as 

( / d^kuj{k)aik, +)a{k, -))^ = {Ho)f = 0(6-^) . (3.6) 

Let {■)t be the time-evolved measure at time t. Its rescaled Wigner function is 

W%y, k, t) = e^ [ dr]e'^^y\a{k - £7^/2, +)a{k + er^/2, -))t/e • (3.7) 

J{T/£)3 



Kinetic theory claims that 



YimW'{[r\,,k,t) = W{r,k,t), (3.^ 



where W{r, k,t) is the solution of the phonon Boltzmann equation 
^W{r, k, t) + ^Vuj{k) ■ VrW{r, k, t) 

ot 277 



y / dkidk2{uj{k)uj{ki)u{k2))-^S{uj{k) + (Tiuj{ki) + a2Uj{k2)) 

(Tl,(T2=±i 



bik + aifci + 02k2) {W{r, ki, t)W{r, fca, t) 

+GiW{r, k, t)W{r, k2, t) + a2W{r, k, t)W{r, k^, t)) (3.9) 

to be solved with the initial condition W{r, k, 0) = W^{r, k). 

The free streaming part is an immediate consequence of the evolution of W as 
generated by Hq. The strength of the cubic nonlinearity was assumed to be of order 
Y^, which results in an effect of order 1 on the kinetic time scale. The specific form 
of the collision operator will be explained in the following section. It can be brought 
into a more familiar form by performing the sum over 0"i,(J2. Then the collision 
operator has two terms. The first one describes the merging of two phonons with 
wave number k and ki into a phonon with wave number k2 = k + ki, while the second 
term describes the splitting of a phonon with wave number k into two phonons with 
wave numbers ki and k2, k = ki + k2. In such a collision process energy is conserved 
and wave number is conserved modulo an integer vector. 

In fj4.7|) the summand with cti = 1 = (T2 vanishes trivially. However it could be 
the case that the condition for energy conservation, 

uj{k) + uj{k') -uj{k + k') = 0, (3.10) 

has also no solution. If so, the collision operator vanishes. In fact, for nearest 
neighbor coupling only, a(0) = 6, a{e) = —1 for |e| = 1, a{x) = otherwise, it can 
be shown that ()3.10p has no solution whenever cuq > 0. To have a non-zero collision 
term we have to require 

dk [ dk'6{u{k) + uo{k') - u{k + k')) > , (3.11) 



which is an implicit condition on the couplings a{x). A general condition to ensure 
()3.11|) is not known. A simple example where ()3.11|) can be checked by hand is 

3 

u;(A;) = cjo + J](l - cos(277A;")) , k = {k\k^,k^). (3.12) 

Q = l 

It corresponds to suitable nearest and next nearest neighbor couplings. 
There is a second more technical condition which requires that 

sup / dk'6{uj{k) + uj{k') - uj{k + k')) = cq < 00 . (3.13) 

k J 

It holds for the dispersion relation (j3.12p . This uniform bound allows for a simple 
proof that the Boltzmann equation has a unique solution for short times provided 
W^{r, k) is bounded. 
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4 Feynman diagrams 

Denoting by the average with respect to the measure at time t (in microscopic 
units), the starting point of the time-dependent perturbation series is the identity 

n n n 

{Y[a{kj,(Tj))t = exp [it{'^(TjUj{kj))](Y[a{kj,aj))^ 

j=i j=i j=i 

+i\/e I ds exp [i(t — s) ( ajUj{kj))~\ 
Jo 

n 

(E E " 

1=1 a',a"=±l 
n 

(( n a{k,,a,))a{k\a')a{k'\a"))^. (4.1) 

Here 

(P{k, k\ k") = X{8uj{k)uj{k')uj{k"))~'^^ (4.2) 

One starts with n = 2 and ai = 1, a2 = —1. Then on the right hand side of (j4.1|) 
there is the product of three a's. One resubstitutes ()4.H) with n = 3, etc. Thereby 
one generates an infinite series, in which only the average over the initial Gaussian 
measure (■)^ appears. 

To keep the presentation transparent, let me assume that (■)^ is a translation 
invariant Gaussian measure with 

{a{k, ±))^ = , {a{k, a)a{k' , a))^ = , 

{a{k, +)a{k', -))^ = 6{k - k')W{k) . (4.3) 

Then the measure at time t is again translation invariant. Kinetic scaling now merely 
amounts to considering the long times t/e. The Wigner function at that time is then 
represented by the infinite series 

oo 

{a{q, -)a{p, = 6{q - p) (w{q) + ^ W^^iq, t)) . (4.4) 

n=l 

The infinite sum is only formal. Taking naively the absolute value at iteration 2n 
one finds that 

\W^iq, t)\ < £"(t/£)2"((2n)!)-i(2n)!c2"((2n + 2)\/2^^\n + 1)!) . (4.5) 

Here = (y^)^", {t / e)"^^ / {2n)\ comes from the time integration, {2n)\ from the 
sum over i in ()4.H) . c^" from the fc-integrations and the initial W{k), and the factor 
{2n + 2)!/2"+^(n + 1)! from the Gaussian pairings in the initial measure. Thus even 
at fixed e there are too many terms in the sum. 



/ dk'dk"(j){ke, k', k")5{-aih + cr'k' + a"k") 
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Since no better estimate is available at present, we concentrate on the structure 
of a single summand W^{q, t). 5{q—p)W^{q^ t) is a sum of integrals. The summation 
comes from 

- the sum over a', a" in ()4.1|) 
-the sum over £ in (j4.H) 

- the sum over all pairings resulting from the average with respect to the initial 
Gaussian measure (■)*^. 

Since each single integral has a rather complicated structure, it is convenient to 
visualize them as Feynman diagrams. 

A Feynman diagram is a graph with labels. Let us first explain the graph. The 
graph consists of two binary trees. It is convenient to draw them on a "backbone" 
consisting of 2n + 2 equidistant horizontal level lines which are labelled from 
(bottom) to 2?2 + 1 (top). The two roots of the tree are two vertical bonds from 
line 2n + 1 to level line 2n. At level m there is exactly one branch point with two 
branches in either tree. Thus there are exactly 2n branch points. At level there 
are then 2n + 2 branches. They are connected according to the pairing rule, see 
Figure below. 

In the Feynman graph each bond is oriented with arrows pointing either up 
(cr = +1) or down (a = —1). The left root is down while the right root is up. If 
there is no branching the orientation is inherited from the upper level. At a pairing 
the orientation must be maintained. Thus at level a branch with an up arrow can 
be paired only with a branch with a down arrow, see ()4.3|) . Every internal line in 
the graph must terminate at either end by a branch point. Every such internal line 
admits precisely two orientations. 
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Next we insert the labels. The level lines to 2ra + 1 are labelled by times 
< ti . . . < < t. The left root carries the label q while the right root carries the 
label p. Each internal line is labelled with a wave number k. 

To each Feynman diagram one associates an integral through the following steps. 

(i) The time integration is over the simplex < ti . . . < t2n. < ^ as dti . . . dt2n- 

(ii) The wave number integration is over all internal lines as j dki . . . J dk^, where 
K = Sra — 1 is the number of internal lines. 

(iii) One sums over all orientations of the internal lines. 
The integrand is a product of three factors. 

(iv) There is a product over all branch points. At each branchpoint there is a root, 
say wave vector ki and orientation ai, and there are two branches, say wave vectors 
^2, ^3 and orientations ct2, (t^. Then each branch point carries the weight 

6{-criki + 0-2^:2 + (r3k3)ai(t){ki, /c2, ^3) • (4.6) 

If one regards the wave vector A; as a current with orientation a, then ()4.6|1 expresses 
Kirchhoff's rule for conservation of the current. 

(v) By construction each bond carries a time difference tm+i ~ t^, a wave vector k, 
and an orientation cr. Then to this bond one associates the phase factor 

exp[i{tm+i - tm)(Tujik)/e] . (4.7) 

The second factor is the product of such phase factors over all bonds. 

(vi) The third factor of the integrand is given by 

n+l 

l[W{k,), (4.8) 

j=i 

where ki, . . . , kn+i are the wave numbers of the bonds between level and level 1. 

(vii) Finally there is the prefactor (—!)"£:"". 

To illustrate these rules we give an example for n = 2, see Figure above. The 
associated integral is given by, more transparently keeping the 5-functions from the 
pairings, 

/ dti . . . dt^ / dki . . . dks 

Jo<ti<...<t4<t JT24 

5{q + ki - k2)6{k2 + k3 + k4)5{-p - k^ - kQ)5{-ki + kj - ks) 
(t){q, ki, ^2)0(^2, ks, kA)(t){p, kr,, ke)(l){ki, kr, ks) 
5{k, - k,)W{k,)5{ks - k3)W{ks)5{k^ - k^)W{k^) 
exp [{i{t - ti){-uOq + uop) + i{ti - t3){uoi - UJ2 + ujp) 

+i{t3 - t2){uJi +003 + 004 + Up) + i{t2 - ti){ljJi + CJ3 + CJ4 - Go's - UJq) 

+itl{uJ^ - + cjg + - - uj&)]/e\ (4.9) 
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with ujq = u{q), Up = uj(p), toj = uj{kj). 

S{q — p)W^{q,t) is the sum over all Feynman diagrams with 2n + 2 levels and 
thus is a sum of oscillatory integrals. In the limit e ^ only a few leading terms 
survive while all remainders vanish. E.g., the Feynman diagram above is subleading. 
In fact, the conjecture of kinetic theory can be stated rather concisely: 

Kinetic Conjecture: In a leading Feynman diagram the Kirchhoff rule never 
forces an identification of the form 6{kj) with some wave vector kj. In addition, the 
sum of the 2{n — m + 1) phases from the bonds between level lines 2m and 2m + 1 
vanishes for every choice of internal wave numbers. This cancellation must hold for 
m = 0, . . . ,n. 

Since we assumed the initial data to be spatially homogeneous, the phonon Boltz- 
mann equation simplifies to 

d 



dkidk2(j){k, ki, /C2) S{u{k) + aiuj{ki) + a2io{k2)) 

T6 



S{k + aiki + 0-2^2) {W{ki, t)W{k2, t) + 2a2W{k, t)W{ki, t)) , (4.10) 

where we used the symmetry with respect to (/ci,(Ji) and (/i;2,o"2)- To ()4.10|) we 
associate the Boltzmann hierarchy 

d 

Q^fn = Cn,n+lfn+l , U = 1,2, . . . , (4.11) 

acting on the symmetric functions fniki, . . . , kn) with 

n „ 

Cn,n+ifn+i{ki,. . . ,kn) = A-kX^^^ ^ / dk' dk" (j){kl, k' , k"Y 

5{uj{h) + a'uj{k') + a"uj{k"))5{h + (r'k' + a"k") 

[fn+i{k,, ...,k', h+i, k") + 2a"f{k,, ...,k^, k')] . (4.12) 

Under the condition (jH.lHj) and provided ||W^||oo < oo, the hierarchy ()4.1H1 has a 
unique solution for short times. In case 

n 

fn{ki,...,kn,0) = l[W{kj), (4.13) 

i=i 

the factorizatin is maintained in time and each factor agrees with the solution of 
the Boltzmann equation ()4.1()j) . From ()4.1H1 one easily constructs the perturbative 
solution to ()4.10|) with the result 

oo _ 

W{k, t) = W{k) + V -r(Ci,2 . . . Cn^n+lW^""^^) {k) 

^ — ^ n\ 

n=l 
oo 

= Ty(A;) + J]iy„(A;,t). (4.14) 



n=l 
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The series in (14 .1411 converges for t sufficiently small. 

For n = 1,2 the oscillating integrals can be handled by direct inspection with 
the expected results lim^^Q Wf{k,t) = Wi{k,t), lim^^o W^2^(^5 ^) = W2{k,t). If the 
leading terms are as claimed in the Kinetic Conjecture, then they agree with the 
series ()4.14|) . The complete argument is a somewhat tricky counting of diagrams, 
which would lead us too far astray. Thus the most immediate project is to establish 
that all subleading diagrams vanish in the limit e — > 0. This would be a step further 
when compared to the investigations 05 0- 

Of course a complete proof must deal with the uniform convergence of the series 
in (lOjl . 

Acknowledgements. I thank Jani Lukkarinen for instructive discussions and Gi- 
anluca Panati for a careful reading. 
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